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ULTRA-HIGH-RESOLUTION OBSERVATIONS 
OF MHD WAVES IN PHOTOSPHERIC 
MAGNETIC STRUCTURES 
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Here we review the recent progress made in the de¬ 
tection, examination, characterisation and interpre¬ 
tation of oscillations manifesting in small-scale mag¬ 
netic elements in the solar photosphere. For white- 
light observations, the solar photosphere is defined 
as the layer that emits most of the Sun’s electro¬ 
magnetic radiation, and therefore creates the limit¬ 
ing depth of optical imaging sequences. This region 
of the Sun’s atmosphere is especially dynamic, con¬ 
stantly evolving, and importantly, permeated with an 
abundance of large- and small-scale concentrations of 
magnetic fields. Such magnetic features can span di¬ 
ameters of hundreds to many tens of thousands of 
km, and are thus commonly referred to as the ‘build¬ 
ing blocks’ of the magnetic solar atmosphere. How¬ 
ever, it is the smallest, and by far the most numerous 
magnetic elements that have risen to the forefront of 
solar physics research in recent years. Structures, 
which include magnetic bright points, are often at 
the diffraction limit of even the largest of solar tele¬ 
scopes. Importantly, it is the improvements in facil¬ 
ities, instrumentation, imaging techniques and pro¬ 


cessing algorithms during recent years that have al¬ 
lowed researchers to examine the motions, dynamics 
and evolution of such features on the smallest spa¬ 
tial and temporal scales to date. It is clear that 
while these structures may demonstrate significant 
magnetic field strengths, their small sizes make them 
prone to the buffeting supplied by the ubiquitous sur¬ 
rounding convective plasma motions. Here, it is be¬ 
lieved that magnetohydrodynamic waves can be in¬ 
duced, which propagate along the field lines, carrying 
energy upwards to the outermost extremities of the 
solar corona. Such wave phenomena can exist in a 
variety of guises, including fast and slow magneto¬ 
acoustic modes, in addition to Alfven waves. Cou¬ 
pled with rapid advancements in magnetohydrody¬ 
namic wave theory, we are now in an ideal position to 
thoroughly investigate how wave motion is generated 
in the solar photosphere, which oscillatory modes are 
most prevalent, and the role that these waves play in 
supplying energy to various layers of the solar atmo¬ 
sphere. 


1. INTRODUCTION 

Ever since oscillatory motion was first discovered 
in the Sun’s tenuous atmosphere [Leighton, 1960; 
Leighton et ah, 1962; Noyes Sz Leighton, 1963], it 
has been a goal amongst physicists to detect, iden¬ 
tify, characterise and understand the diverse va¬ 
riety of wave modes manifesting in the solar at¬ 
mosphere. Initial observations in the optical por¬ 
tion of the electromagnetic spectrum allowed os¬ 
cillations present in the photospheric and chro¬ 
mospheric layers to be studied. However, during 
these early stages technology was not as advanced 


Astrophysics Research Centre, School of 
Mathematics and Physics, Queen’s University 
Belfast, Belfast, Northern Ireland, BT7 INN, UK. 

2 Solar Physics and Space Plasma Research 
Centre (SP 2 RC), The University of Sheffield, Hicks 
Building, Hounsfield Road, Sheffield, S3 7RH, UK. 


as it is today, and as a result, key modern tech¬ 
niques such as adaptive optics [AO; e.g., Rinnnele 
&; Marino, 2011], multi-object multi-frame blind 
deconvolution [MOMFBD; van Noort et ah, 2005] 
and speckle reconstruction [Woger et ah, 2008] 
were unavailable to help combat the fine-scale im¬ 
age degradation caused by the Earth’s atmosphere. 
Therefore, initial research was dedicated to probing 
large-scale solar structures, including sunspots and 
super-granules. However, even at these large spa¬ 
tial scales, a wealth of oscillatory phenomena were 
found to be omnipresent [e.g., Deubner & Liedler, 
1969; Ulrich, 1970; Deubner, 1971]. Such waves 
demonstrated periodic intensity and velocity fluc¬ 
tuations, and thus were placed under the same um¬ 
brella as acoustic modes, which are ultimately de¬ 
fined by their intrinsic signatures of compressions 
and rarefactions. Difficulties were encountered in 
follow-up work when the measured phase velocities 
of the waves were found to be too large to be ex¬ 
plained by over-simplistic acoustic modelling. In- 
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stead, Osterbrock [1961] and Mein & Mein [1976] 
hypothesised that the magnetic fields in which the 
waves were embedded within must also be taken 
into consideration. 

Over the following decades the examination of 
magneto-acoustic waves rose to the forefront of ob¬ 
servational solar physics, with Ulmschneider [1976] 
aptly asking whether they may provide a signif¬ 
icant channel for energy leaking into the outer 
layers of the Sun’s atmosphere. Furthermore, 
the mere presence of a magnetic field also intro¬ 
duces a number of additional viable wave modes 
that may simultaneously exist within the atmo¬ 
sphere. Through magneto-hydrodynamic (MHD) 
modelling it was found that although these addi¬ 
tional waves have similarities with purely acoustic 
modes, they are often highly anisotropic. This is 
because the addition of the magnetic field intro¬ 
duces a dependency on both ( i ) the alignment of 
the wavevector, k, with the direction of the back¬ 
ground magnetic field, Bq, and ( ii ) the ratio of 
the kinetic pressure, po, to the magnetic pressure, 
Hg/2/io (or Bq/8it in cgs units), in the environ¬ 
ment which supports the wave. This ratio is com¬ 
monly referred to as the plasma /3, defined as (3 = 
2poPo/L>o, where po is the magnetic permeability 
of free space. This quantity can be re-written in 
terms of the local Hydrogen number density, njj, 
the plasma temperature, T, and the Boltzmann 
constant, ks, as /3 = SimuTkE/B^, providing a 
representation of the plasma (5 in cgs units. Thus, 
in the lower layers of the solar atmosphere where 
the temperature is relatively low (T ~ 6000 K) 
and the magnetic field strength is still intensely 
concentrated (Bq > 1000 G), it is clear to see how 
a majority of the magnetic structures supporting 
magneto-hydrodynamic wave phenomena are of¬ 
ten defined by /3<Cl (i.e., dominated by magnetic 
pressure). This has important consequences for 
the wave modes that are expected to exist within 
this plasma regime, and allows for the manifesta¬ 
tion of ‘fast’ and ‘slow’ magneto-acoustic waves, 
in addition to Alfven waves [Edwin & Roberts, 
1983]. Structures in the lower solar atmosphere 
which demonstrate oscillatory behaviour are often 
observed to be long, elongated features (e.g., mot¬ 
tles, spicules, fibrils, etc.), and as such are typically 
modelled by employing cylindrical geometry. Here, 
an overdense magnetic ‘flux tube’ will create an ef¬ 
ficient waveguide that will support an even richer 
variety of MHD waves depending on the azimuthal 
wavenumber, m, including sausage (rn = 0), kink 
(m = 1) and fluting modes (m >1), which are yet 
further classified via their trapped/leaky, fast/slow 
and body/surface characteristics. Thus, the quest 
to detect all of the various MHD wave modes has 
intensified, but so too has the drive to determine 


the precise role they play in the transportation 
and dissipation of energy through the Sun’s atmo¬ 
sphere. 

The ongoing research into waves and oscillations 
in the solar atmosphere does not take place for 
solely esoteric purposes. Rather, it has the ability 
to allow us to delve deep below the visible solar sur¬ 
face through helioseismology techniques [e.g., Du¬ 
vall et al., 1993; Schou et al., 1998; Lopes & Silk, 
2014], and to understand the coupling between 
(quasi-) periodic flows and plasma motions that 
are abundantly apparent over a wide range of spa¬ 
tial scales (e.g., spicules, mottles, fibrils, plumes, 
prominences, etc.). Quantifying the possible con¬ 
tribution of observed MHD waves to plasma heat¬ 
ing is also of immense interest in solar physics. 
The paradoxical nature of how the Sun’s outer at¬ 
mosphere is heated to (and maintained at) multi¬ 
million degree temperatures is a problem that has 
been plaguing scientists for over half a century. 
Even the solar chromosphere, a thin atmospheric 
layer that is only heated to a few thousand de¬ 
grees above the underlying photosphere, requires 
extraordinary plasma heating processes to balance 
the radiative losses experienced in this relatively 
high-density environment. It has long been be¬ 
lieved that MHD waves, generated in the pho¬ 
tosphere and channelled upwards along magnetic 
field lines, may be able to contribute directly to 
plasma heating providing a suitable (and efficient) 
conversion mechanism exists. 

To generate localised heating from MHD waves, 
smaller length scales must be created via physical 
processes such as resonant absorption, phase mix¬ 
ing or shock formation. Observing such fine-scale 
plasma dynamics directly presents a significant ob¬ 
servational challenge due to the near- (or sub-) res¬ 
olution scales involved. De Pontieu et al. [2004] 
demonstrated how ubiquitous p-rnode oscillations 
can be channelled along magnetic field lines, ulti¬ 
mately giving rise to dynamic phenomena in the 
chromosphere and transition region. Furthermore, 
Hansteen et al. [2006] revealed how shocks can form 
when slow magneto-acoustic waves, which are gen¬ 
erated by convective flows and global p-rnode oscil¬ 
lations, leak upwards from the solar surface along 
magnetic field lines and encounter the steep den¬ 
sity gradients intrinsic to the chromosphere. Of 
significant importance is the fact that there is now 
overwhelming evidence to suggest that waves have 
the ability to deform magnetic field lines and in¬ 
duce the necessary instabilities required to incite 
reconnective phenomena over a wide range of at¬ 
mospheric heights [e.g., Isobe &; Tripathi, 2006; 
Isobe et al., 2007; Jess et al., 2010a; Li Sz Zhang, 
2012; Jackiewicz & Balasubramaniam, 2013; Shen 
et al., 2014]. Thus, examining the generation and 
behaviour of oscillatory phenomena in the solar 



D. B. JESS & G. VERTH: MHD WAVES IN THE SOLAR PHOTOSPHERE • 3 



Figure 1. A G-band image (left), acquired by the ROSA multi-camera system at the Dunn Solar Tele¬ 
scope, revealing a large number of MBPs. These features are visible as intensity enhancements within 
the intergranular lanes, and often demonstrate magnetic field strengths exceeding 1000 G. The axes are 
in heliocentric arcseconds, where 1" ~ 0.725 Mm. The red square highlights a 12 x 12 Mm 2 subregion 
of the field-of-view, with an equally-sized simulated G-band image displayed in the upper-right panel. 
The output of the MuRAM numerical code also reveals bright magnetic features manifesting within 
the intergranular lanes, suggesting a significant amount of agreement between current high-resolution 
observations and radiative magnetohydrodynamic modelling techniques. The lower-right panel displays 
the local sound speeds calculated from the local densities and pressures in the simulated polytropic at¬ 
mosphere. It is clear that intergranular lanes, where magnetic bright points reside, often display sound 
speeds on the order of 10 kms . Images adapted from Jess et al. [2012c]. 


atmosphere has the potential to shine light on a 
wide variety of physical phenomena, from magnetic 
reconnection on sub-arcsecond scales through to 
prominence eruptions spanning many hundreds of 
arcseconds. 

In more modern times, the advent of high- 
sensitivity and low-noise camera systems has al¬ 
lowed short-exposure and high-cadence image se¬ 
quences to be obtained from both ground- and 
space-based observatories. In particular, for 
ground-based solar telescopes such as the Dunn 
Solar Telescope (DST) and the Swedish Solar Tele¬ 
scope (SST), the successful application of adaptive 
optics and post-processing techniques has allowed 
researchers to focus their attention on the small¬ 
est magnetic elements observed in the Sun’s lower 
atmosphere. Features that are close to the res¬ 
olution limits of current solar telescopes, includ¬ 
ing magnetic bright points [MBPs; Dunn Sz Zirker, 


1973; Stenflo, 1985; Solanki, 1993], offer unique 
advantages over more large-scale structures (e.g., 
sunspots) in the study of MHD wave phenom¬ 
ena. For instance, these omnipresent elements, eas¬ 
ily identifiable and ideal for feature tracking algo¬ 
rithms, are not as magnetically complex as macro¬ 
sized structures such as sunspots, and are there¬ 
fore more readily compared with cylindrical geom¬ 
etry MHD approximations. Furthermore, due to 
their small size [often less than 0.3 ,/ or 220 km in 
diameter; Crockett et al., 2010], these structures 
are more prone to the buffeting imposed by the 
complex evolution of surrounding granules, thus 
increasing the likelihood of generating wave mo¬ 
tion at the photospheric base of the magnetic field 
lines. In addition, the magnetic field strengths as¬ 
sociated with MBP structures are of comparable 
magnitude to those found in large-scale sunspots 
[> 1 kG; Cauzzi et al., 2000; Jess et al., 2010b], 
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making them extremely viable conduits for carry¬ 
ing a variety of MHD waves, including fast and slow 
magneto-acoustic modes, in addition to Alfven 
waves. MBPs have demonstrated their ubiquity 
in both observational and simulated datasets, with 
an example of such features displayed in Figure 1. 

The simultaneous development of theoretical, 
numerical and analytical modelling tools for the 
lower solar atmosphere has allowed the unifica¬ 
tion of MHD with seismological techniques [e.g., 
Kim et ah, 2008; Verth et al., 2011; Kuridze et 
ah, 2013; Morton et ah, 2014], Now, observers 
are able to quantify and understand oscillatory pa¬ 
rameters that may be below the resolution limit 
imposed by even the largest of modern ground- 
based optical telescopes (e.g., the 1.6m New Solar 
Telescope, NST, at the Big Bear Solar Observa¬ 
tory). This has important consequences, particu¬ 
larly when attempting to identify multiple modes 
existing within the same structure [e.g., Morton 
et ah, 2012], or when trying to diagnose small- 
amplitude waves that may often become swamped 
by instrumental noise and/or the point spread 
function of the telescope. In this chapter, we will 
review the recent observations of waves and oscilla¬ 
tions manifesting in fine-scale magnetic structures 
in the solar photosphere, which are often inter¬ 
preted as the “building blocks” of the magnetic 
Sun. 

2. MAGNETO-ACOUSTIC WAVES 

The launch of the Hinode [Kosugi et ah, 2007] 
space telescope, equipped with the 0.5 m Solar 
Optical Telescope [SOT; Suematsu et ah, 2008; 
Tsuneta et ah, 2008], provided researchers with 
easily accessible and high resolution (~0.2 ,/ or 
150 km) observations of the lower solar atmo¬ 
sphere. The lack of atmospheric turbulence results 
in the images being ‘seeing free’, and thus suit¬ 
able for long duration studies of small-scale mag¬ 
netic elements on the surface of the Sun. Notably, 
Carlsson et ah [2007] were at the forefront of em¬ 
ploying the high resolution Ca II H and contin¬ 
uum filters on-board SOT to examine the propa¬ 
gation of acoustic waves down to the diffraction 
limit of the instrumentation. Through examina¬ 
tion of all spatial locations (both magnetic and 
non-magnetic) Carlsson et ah [2007] speculated 
that the total energy flux provided, even at such 
high spatial and frequency resolutions, was insuffi¬ 
cient to contribute to atmospheric heating. How¬ 
ever, Wedemeyer-Bohm et ah [2007] stated that 
the methods used may overlook dynamic patterns 
created on sub-resolution scales, and as a result 
severely underestimate the actual mechanical flux 
[Kalkofen, 2007, 2008]. Importantly, this work 


inspired many other solar scientists to focus on 
highly magnetic photospheric elements, where the 
strong magnetic flux concentrations may promote 
more efficient energy propagation. 

Employing high spatial and temporal resolution 
observations from the ground-based DST, Jess et 
ah [2007] detected more prevalent photospheric os¬ 
cillations in what was believed to be magnetic con¬ 
centrations surrounding a large sunspot. These os¬ 
cillations demonstrated power well in-excess of the 
quiescent background, but unfortunately no direct 
information on the associated magnetic fields were 
available, and thus the waves were tentatively tied 
to magneto-acoustic phenomena. Further work by 
Andie [2007] found that the locations of small- 
scale oscillatory power in the photosphere corre¬ 
lated well with red-shifted velocities (i.e., down¬ 
flows). This is in agreement with the theory that 
many small-scale magnetic elements in the pho¬ 
tosphere are formed via the process of convective 
collapse [Spruit, 1976], which is further substan¬ 
tiated by the modern-day observations and MHD 
simulations of MBPs [e.g., Utz et ah, 2013; Hewitt 
et ah, 2014]. However, it must be stated that while 
such oscillations have been tied to MBPs display¬ 
ing red-shifted (i.e., downflow) velocities, this does 
not necessarily mean that the embedded waves are 
also downwardly propagating. A bulk motion may 
exist within the confines of the magnetic flux tube, 
hence giving rise to red-shifted Doppler signatures, 
yet the phase velocity of the wave may have a suffi¬ 
ciently large upward magnitude such that the over¬ 
all group velocity, which crucially describes the ve¬ 
locity at which energy is propagated by the wave, is 
directed upwards. For example, the observational 
work of Andie [2007] and Narayan [2011], corrobo¬ 
rated by the numerical modelling presented by He¬ 
witt et al. [2014], established downflow velocities 
on the order of a few kms -1 within the confines 
of small-scale photospheric magnetic elements (see, 
e.g., the upper panels of Figure 2). These subsonic 
plasma flows may easily become overshadowed by 
upwardly propagating magneto-acoustic phase ve¬ 
locities that are close to the photospheric sound 
speed [~10 kms -1 ; Jess et al., 2012c]. As a result, 
the group velocity, and therefore the direction of 
energy propagation, would be directed upwards. 
Recently, Kato et al. [2011] further revealed how 
the processes intrinsic to convective collapse may 
also drive magneto-acoustic wave phenomena in 
small-scale magnetic elements, irrespective of the 
presence of sub-surface p-mode oscillations. The 
authors employed radiative MHD simulations and 
demonstrated how the coupling between external 
downdrafts in the intergranular lanes and the mo¬ 
tions of the embedded plasma relies heavily on the 
inertial forces that act on the magnetic flux con¬ 
centration. These forces act to ‘pump’ the internal 
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Figure 2. Simultaneous G-band intensity (upper left) and line-of-sight velocity (upper right) sub-fields 
extracted from a synthesised field-of-view created using the MuRAM radiative magnetohydrodynamic 
code. Each image covers less than 1 square arcsecond, and is centred on a recently formed MBP structure 
manifesting within the intergranular lanes of the surrounding convective plasma. The red contours in 
the upper-left panel outline locations where the downflow (i.e., red-shifted) velocities exceed 3 kms” , as 
defined by the graduated colour spectrum displayed at the right-hand side of the Figure. The red-shifted 
plasma encompasses both the MBP and the surrounding intergranular lanes, indicating the process of 
convective collapse is likely to play an important role in the formation of MBP features. Note that the 
downflow velocities at the centre of the MBP, where wave phenomena is likely to manifest as a result of 
the increased magnetic field strengths, are slightly weaker than at its perimeter. The lower panels display 
the occurrence of oscillations simultaneously visible in ROSA G-band and 4170 A continuum images, as 
a function of the oscillation period and phase angle for regions containing MBPs (lower left) and those 
without (lower right). The colour scale represents the number of detections as a percentage of the total 
events with an associated coherence level exceeding 85%. A horizontal dashed line represents a phase 
angle of 0 degrees, while dotted lines highlight a region inside which detections become unreliable due 
to cadence restrictions (0.528 s for the observational time series). A preference for negative phase shifts 
highlights the abundance of upwardly propagating wave motion in non-magnetic and magnetic photo- 
spheric features, including the small-scale MBPs shown in the upper panels. Upper images courtesy of 
R.L. Hewitt, and based on the data presented by Hewitt et al. [2014], while the lower panels have been 
adapted from Jess et al. [2012c], 


atmosphere of the magnetic flux tube in a down¬ 
wards direction, which eventually causes the at¬ 
mosphere to rebound, producing upwardly propa¬ 
gating magneto-acoustic waves along the magnetic 
field lines [Kato et al., 2011]. 

Through examination of the magneto-acoustic 
wave dynamics associated with large-scale sunspots, 
Nagashima et al. [2007] employed Hinode/SOT ob¬ 
servations to reveal how oscillatory power is often 
drastically reduced in the presence of strong mag¬ 
netic field concentrations; a common phenomenon 


now referred to as ‘acoustic power suppression’ 
[e.g., Woods & Cram, 1981; Thomas et al., 1982; 
Title et al., 1992; Parchevsky & Kosovichev, 2007; 
Chou et al., 2009; Ilonidis &: Zhao, 2011; Cou- 
vidat, 2013]. Lawrence & Cadavid [2010, 2012] 
and Chitta et al. [2012a] were able to corrobo¬ 
rate these general findings through examination 
of Hinode/SOT G-band image sequences contain¬ 
ing MBPs. However, interestingly Chitta et al. 
[2012a] found evidence to suggest that magneto¬ 
acoustic power at the highest temporal frequen- 
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cies (i.e., periodicities less than 100 s) actu¬ 
ally demonstrated power amplification. Follow-up 
work, employing the Rapid Oscillations in the So¬ 
lar Atmosphere [ROSA; Jess et ah, 2010c] multi¬ 
camera imaging system on the DST, revealed how 
magneto-acoustic waves with periodicities below 
100 s demonstrate significant turbulent compo¬ 
nents within their power spectra [Lawrence et ah, 
2011]. Lawrence et ah [2011] suggest that the ob¬ 
served magneto-acoustic waves may be generated 
by the interaction of the magnetic held lines with 
plasma downflows (i.e., characteristic of the con¬ 
vective collapse process) that are very turbulent in 
their nature. 

While the observed periodicities of magneto¬ 
acoustic waves in MBPs generally span the entire 
p-mode spectrum, there is increasing evidence to 
suggest that the underlying magnitude of the mag¬ 
netic held directly influences the dominant period. 
Kostik & Khomenko [2013] employed the Triple 
Etalon SOlar Spectrometer [TESOS; Tritschler et 
ah, 2002] on the German Vacuum Tower Telescope 
(VTT) to obtain high-resolution spectroscopy of 
the photospheric Ba n absorption line. They found 
that the dominant period of oscillations increases 
by 15-20% as the local magnetic held strength in¬ 
creases from 500 to 1500 G. This has important 
implications since it suggests that the strong mag¬ 
netic fluxes inherent to MBPs may be able to as¬ 
sist the propagation of lower frequency (i.e., below 
3 mHz) magneto-acoustic waves into the chromo¬ 
sphere and beyond, especially if the magnetic held 
lines are suitably inclined to reduce the impact of 
the acoustic cut-off frequency [i.e., Bel & Leroy, 
1977]. As suggested by De Pontieu et ah [2004], de 
Wijn et ah [2009] and Stangalini et ah [2011], to 
name but a few, the efficient leakage of lower fre¬ 
quency magneto-acoustic wave modes into the up¬ 
per solar atmosphere may be able to drive a wide 
variety of high-temperature phenomena, including 
the oscillations observed in coronal fans, plumes 
and loops [e.g., Deforest & Gurman, 1998; Ofman 
et ah, 1999; Ofman & Wang, 2002; De Moortel & 
Hood, 2003, 2004], 

High resolution observations have clearly shown 
the existence of magneto-acoustic wave phenom¬ 
ena in small-scale photospheric magnetic elements. 
However, it is also of significant importance to de¬ 
termine whether these oscillations are propagating 
or standing waves from the viewpoint of supply¬ 
ing energy to the upper levels of the solar atmo¬ 
sphere. Jess et ah [2012c] employed the ROSA 
system with blue continuum (4170 A) and G-band 
filters to examine the propagation of waves be¬ 
tween two discreet layers in the lower solar atmo¬ 
sphere. The MuRAM radiative magnetohydrody¬ 
namic code [Vogler et ah, 2005] was utilised to de¬ 


termine the formation heights of the two ROSA 
filtergrams through a comparison of their corre¬ 
sponding response functions. It was found that 
the continuum and G-band images were separated 
by 75 km in height, thus allowing any propaga¬ 
tion of waves between these two discreet layers to 
be investigated through phase-difference analysis. 
Jess et ah [2012c] were able to detect a wealth 
of oscillatory phenomena in both continuum and 
G-band images. However, only oscillations with 
periodicities above ~140s demonstrated coherent 
phase delays between the adjacent bandpasses. 
A —8° phase lag (traversing a physical displace¬ 
ment of ~75 km) indicated upwardly propagating 
phase speeds on the order of 8 kms -1 . These ve¬ 
locities are similar to the expected sound speed 
(lower-right panel of Figure 1), and thus demon¬ 
strated the linear nature of magneto-acoustic wave 
phenomena in the lower solar atmosphere. Fur¬ 
thermore, as revealed in the lower panels of Fig¬ 
ure 2, the authors determined that 76% of all 
MBP structures demonstrated upwardly propagat¬ 
ing magneto-acoustic wave signatures, helping to 
explain why the outer regions of the solar atmo¬ 
sphere are so ubiquitously populated with MHD 
wave phenomena. Indeed, extensive work has re¬ 
cently been implemented to uncover the connectiv¬ 
ity between lower atmospheric propagating waves 
and running oscillations ubiquitously observed in 
the solar corona [e.g., Tomczyk & McIntosh, 2009; 
Jess et ah, 2012a], 

More recently, Andie et al. [2013] employed 
broadband TiO images obtained using the NST 
to examine the connection between photospheric 
oscillations and the dynamic motions of small- 
scale magnetic flux concentrations. The TiO fil¬ 
ter used is centred on an absorption band of 
molecules around 7056.8 A (incorporating a 10 A 
filter width), and thus averages over all inherent 
absorption and continua contributions, causing the 
resulting images to be only weakly dependent on 
the properties of individual spectral lines. Ulti¬ 
mately, this means that the intensity time series 
will be representative of the true solar continuum. 
Furthermore, since the bandpass is approaching 
the near-infrared, the images will be less sensitive 
to atmospheric seeing variability, and thus provide 
better temporal coverage of time series obtained 
during mediocre weather conditions. The authors 
undertook Fourier and Hilbert transformations of 
the TiO time series related to MBPs, in order to 
extract the amplitude and phase relationships of 
the embedded oscillations, and suggested that the 
detected wave motion is likely to be too complex to 
be generated by a single oscillatory source. Nev¬ 
ertheless, Andie et al. [2013] provided direct ev¬ 
idence for the presence of upwardly propagating 
wave trains in the immediate vicinity of red-shifted 
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Figure 3. Panels displaying the time evolution of circularly polarised (i.e., longitudinal) magnetic fields 
in the photosphere, as captured by the IMaX instrument on-board the Sunrise balloon flight. What 
may initially appear to be pore-sized structures are in fact sub-arcsecond magnetic concentrations con¬ 
tained within the intergranular lanes, with white and black colours representing positive and negative 
magnetic flux densities, respectively, computed using the weak-field approximation. The blue contours 
highlight various iso-magnetic flux densities, while the yellow contours represent time-constant magnetic 
fluxes equal to —4.5 x 10 16 Mx (upper row) and —5.0 x 10 16 Mx (lower row). It is clear that the time- 
constant iso-magnetic contours contract and expand with time, suggesting the presence of compressive 
magneto-acoustic waves embedded within the field lines. Images adapted from Martinez Gonzalez et al. 
[ 2011 ]. 


(i.e., downflowing) material, suggesting the phase 
velocities of the magneto-acoustic wave phenom¬ 
ena were significantly higher than 5 kms' 1 at the 
photospheric layer. 

2.1. Sausage Waves 

Even though the sausage mode is the lowest az¬ 
imuthal order compressible mode (i.e., m = 0), it 
has still proved extremely difficult to identify in ob¬ 
servations. Through mathematical understanding 
and numerical modelling, these waves will demon¬ 
strate observational characteristics consistent with 
the simultaneous periodic intensity and area fluctu¬ 
ations of the magnetic flux tube. One of the main 
obstacles, at least observationally when attempting 
to detect sausage-mode oscillations, is a combina¬ 
tion of the instrumental spatial resolution (i.e., to 
be able to detect the fractional area changes) in ad¬ 
dition to the detector sensitivity (i.e., to be able to 
extract the small-scale intensity fluctuations over 
the intrinsic background noise). Only with the con¬ 
sistently high resolving power of modern solar fa¬ 
cilities, coupled with the low noise characteristics 
synonymous with cooled CCD and CMOS detec¬ 
tors, has it been possible to detect sausage waves 
in small-scale photospheric magnetic elements. 

Fujimura &; Tsuneta [2009] revealed the true 
power of the SOT on-board the Hinode spacecraft 
by examining the intensity and velocity oscilla¬ 
tion characteristics in relation to the vector mag¬ 
netic field (see, e.g., the lower panels of Figure 5). 
The authors found, through phase relationships 
between the various waveforms, that small-scale 
MBPs in the photosphere demonstrated signatures 


of specific magneto-acoustic waves, in particular 
the sausage and kink modes. The observed fluc¬ 
tuations in the magnetic field, the line-of-sight ve¬ 
locity and the structure’s intensity indicated root- 
mean-square amplitudes of 4 — 17 G (0.3%—1.2%), 
0.03 — 0.12 kins' 1 , and 0.1%—1%, respectively. 
The small amplitudes of the observed fluctuations 
emphasised the importance of high-resolution (spa¬ 
tial, temporal and spectral) observations when 
attempting to diagnose certain magneto-acoustic 
wave modes. Importantly, the detected oscillations 
maintained significant overlap with the global p- 
mode periodicities, further confirming the assump¬ 
tion that the majority of wave motion found in the 
solar atmosphere is driven by the underlying and 
omnipresent p-rnode oscillations. 

Employing the spectropolarimetric capabilities 
of the Imaging Magnetograph experiment [IMaX; 
Martinez Pillet et al., 2011] oboard the Sunrise bal¬ 
loon flight, Martinez Gonzalez et al. [2011] estab¬ 
lished clear evidence for periodic fluctuations in the 
magnetic field strength associated with small-scale 
(<1 7/ ) elements within the intergranular lanes. The 
isocontours related to specific magnetic fluxes were 
observed to oscillate in both area and magni¬ 
tude, with some antiphase behaviours suggesting 
the presence of sausage-mode waves [Fujimura & 
Tsuneta, 2009]. However, the authors found that 
the periodicities of the oscillations were seldomly 
constant, and in fact often varied by several min¬ 
utes within the same location (Figure 3). As a re¬ 
sult, Martinez Gonzalez et al. [2011] proposed that 
the wave driver may not be the expected p-mode 
oscillations, and instead may be a consequence of 
granular “forcing”. It was proposed that since the 
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average magnetic fields contained within the inter¬ 
granular lanes may have strengths lower than the 
photospheric equipartition field [~300—500 G; Lin, 
1995; Khomenko et al., 2003; Martinez Gonzalez 
et ah, 2008], the continuously buffeting nature of 
granular flows may directly impose periodic field 
amplification and weakening through the processes 
associated with plasma forcing. Importantly, as 
raised by Martinez Gonzalez et al. [2011], is the 
question as to whether-or-not these fluctuations 
(regardless of whether they are driven by under¬ 
lying p-modes or by granular forcing) are able to 
propagate upwards through the solar atmosphere, 
thus facilitating the relocation of energy to higher 
atmospheric heights. Utilising other magnetically 
sensitive absorption lines, particularly those origi¬ 
nating in the chromosphere (e.g., the Ca II infrared 
triplet at 8542 A) , will allow such magnetic fluctu¬ 
ations to be tracked through the solar atmosphere 
and shine new light on whether they can act as an 
efficient energy conduit. 

Inspired by the results of Fujimura & Tsuneta 
[2009], Moreels & Van Doorsselaere [2013a] and 
Moreels et al. [2013b] developed stringent phase 
relationships that allow the characterisation of 
fast/slow, body/surface and standing/propagating 
sausage-mode waves based on the measured de¬ 
lays between the intensity, cross-sectional area and 
velocity components of the plasma. These rela¬ 
tionships were successfully applied to the observa¬ 
tional work of Morton et al. [2011], which demon¬ 
strated their accuracy and suitability for interpret¬ 
ing particular sausage mode properties, e.g., stand¬ 
ing/propagating, fast/slow, body/surface. How¬ 
ever, the work of Morton et al. [2011] examined so¬ 
lar pores, which with diameters >1250 km are sub¬ 
stantially larger than the ~220 km sizes associated 
with the smallest scale magnetic elements in the 
photosphere [Crockett et al., 2010]. Consequently, 
efforts are now being directed towards high reso¬ 
lution observations of MBPs, which will allow the 
phase relationships to be tested more thoroughly 
for the smallest magnetic structures currently re¬ 
solvable. 

Jess et al. [2012b] examined the mode-coupling 
between compressible and incompressible waves 
found in MBPs, and the connection with their 
chromospheric spicule counterparts. The primary 
aim of this work was not to study sausage-mode 
oscillations, but instead their source. Using the 
Lare2D numerical code [Arber et al., 2001] to 
model an MBP as a thin magnetic flux tube, the 
authors found that a 90° out-of-phase behaviour 
of upwardly propagating magneto-acoustic waves 
at the photospheric layer directly incited the gen¬ 
eration of sausage mode oscillations in the same 
flux tube. Jess et al. [2012b] interpreted the nu¬ 
merical output as evidence for how velocity gra¬ 


dients embedded within the flux tube, as a re¬ 
sult of the out-of-phase magneto-acoustic oscilla¬ 
tions, cause the central axis of the magnetic fields 
to displace transversally. In addition, compres¬ 
sions and expansions in the waveguide are simul¬ 
taneously induced, thus promoting the manifes¬ 
tation of both compressible sausage modes and 
incompressible kink waves at upper-photospheric 
heights. The work of Jess et al. [2012b] clearly 
shows how thin, magnetic structures omnipresent 
throughout the solar atmosphere can readily sup¬ 
port sausage-mode wave generation and propaga¬ 
tion, which is in agreement with the chromospheric 
work of Morton et al. [2012], 

2.2. Kink Waves 

The m = 1 kink mode is unique in that it is 
the only value of the azimuthal wavenumber, m, 
that produces a transverse displacement of a mag¬ 
netic flux tube. Hence, it is readily observed in 
over-dense solar structures with imagers of suf¬ 
ficient spatial and/or temporal resolution. The 
kink mode is highly Alfvenic, since its main restor¬ 
ing force is magnetic tension, thereby making it 
only weakly compressible [see, e.g., Goossens et 
al., 2009, for a detailed discussion]. Kink waves 
have been most extensively studied in the corona 
following the launch of the Transition Region and 
Coronal Explorer [TRACE; Handy et al., 1999] 
spacecraft in 1998. Since 2007 they have also be¬ 
come the subject of much interest in the chro¬ 
mosphere due to their heightened visibility in off- 
limb spicules with Hinode. Since propagating kink 
waves in both chromospheric and coronal waveg¬ 
uides are now seen to be ubiquitous, it is widely be¬ 
lieved they are being driven from the wealth of me¬ 
chanical energy permeating the photosphere. As 
a result, kink waves have become hypothesised as 
a favourable transport mechanism for channelling 
energy from the photospheric convective motions 
through to the upper layers of the solar atmosphere 
[e.g., Cranmer &; van Ballegooijen, 2005; Suzuki & 
Inutsuka, 2005; Verdini & Velli, 2007; Matsumoto 
& Shibata, 2010, to name but a few]. Investigations 
into their possible photospheric signatures are now 
gaining momentum. There is some tentative obser¬ 
vational evidence of the kink instability in penum- 
bral filaments surrounding sunspots [e.g., Ryutova 
et al., 2008; Bharti et al., 2012], but thus far there 
has been no conclusive statistical studies under¬ 
taken of kink waves propagating along such struc¬ 
tures. However, since penumbral filaments are pre¬ 
dominantly highly-inclined and/or closed magnetic 
structures confined to the lower atmosphere they 
therefore do not represent the conduits required 
to transfer kink wave energy into the corona. A 
much more promising form of investigation is to 
track the photospheric motions that could be pos- 
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Figure 4. A broadband Ha image (left), acquired by the SST on 2006 June 18, revealing a large assort¬ 
ment of MBPs within the intergranular lanes. Following correlation tracking on each of the 97 detected 
MBPs within the held-of-view, the resulting power spectrum of the horizontal motions is displayed as a 
function of frequency (solid black line) in the right-hand panel. The dashed line represents a standard¬ 
ised Lorentz profile using identical free parameters (e.g., the correlation time) to that measured in the 
SST observations. It is clear that for frequencies >0.02 Hz (<50 s) the observational horizontal motions 
have more power, highlighting the fact that dynamics on short timescales may be very important in the 
generation of energetic kink and/or Alfven wave phenomena. Images adapted from Chitta et al. [2012b]. 


sible drivers which instigate kink waves. Since fine- 
scale chromospheric magnetic structures, such as 
spicules, fibrils and mottles, are rooted in inter¬ 
granular lanes, their footpoints often reveal them¬ 
selves in photospheric G-band filtergrams as MBPs 
[Jess et al., 2012b]. Tracking the horizontal veloc¬ 
ity components of MBPs provides a useful proxy 
for detecting the transverse wave drivers that can 
excite kink modes, assuming that the inclinations 
of the associated magnetic flux tubes are not too 
far from the vertical at photospheric heights. The 
advent of high temporal and spatial resolutions 
from modern facilities has allowed observers to 
measure such small-scale transverse exertions, and 
as a result, the quest is now on to relate such pho¬ 
tospheric horizontal velocity power spectra with 
that of kink oscillations observed higher up in the 
chromosphere and corona. A key scientific goal is 
now to provide an all-encompassing understanding 
of kink wave excitation, propagation and damping 
throughout the entire solar atmosphere. 

In a series of recent observational papers, Stan- 
galini et al. [2013a, b, 2014] employed data, ac¬ 
quired by the Hinode/SOT and the Sunrise/IMaX 
instruments, to examine the interactions between 
MBPs and their surrounding granular environ¬ 
ment in an attempt to uncover how transverse 
waves in the photosphere are generated. Firstly, 
Stangalini et al. [2013a] found an abundance 
of photospheric magnetic elements demonstrating 
buffeting-induced transverse oscillations with pe¬ 


riodicities under 100 s and velocity amplitudes of 
the order of 1—2 kms^ 1 . The authors interpreted 
their results as evidence for kink waves being gen¬ 
erated by granular buffeting, which are accompa¬ 
nied by longitudinal magneto-acoustic oscillations 
being generated via non-linear interactions. The 
results put forward by these authors corroborates 
numerous theoretical and numerical studies that 
suggest that such motions are a natural response of 
the magnetic field to external plasma forcing [e.g., 
Roberts, 1983; Steiner et al., 1998; Hasan et al., 
2003; Musielak & Ulmschneider, 2003; Khomenko 
et al., 2008; Fedun et al., 2011a; Morton et al., 
2014, to name but a few]. Of particular note 
is the related observational work of Keys et al. 
[2011], who found that, at any one time, approxi¬ 
mately 6% of photospheric MBPs display transver¬ 
sal velocities exceeding 2 kms -1 , which according 
to the mathematical analyses of Choudhuri et al. 
[1993a, b] is more than sufficient to effectively and 
efficiently drive kink-mode oscillations in photo¬ 
spheric magnetic elements. 

Chitta et al. [2012b] employed broadband Ha 
images, acquired by the SST with a cadence of 5 s, 
to calculate the flow velocities of MBPs, relative to 
their local environment, through correlation track¬ 
ing algorithms. The broadband nature of the Ha 
filter (8 A full-width half-maximum) resulted in the 
images being dominated by photospheric contri¬ 
butions, as can be seen in the left panel of Fig¬ 
ure 4. The authors were able to detect, track and 
analyse 97 individual MBP features, which allowed 
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Figure 5. A small sub-field of a ROSA G-band image (upper left), obtained using the DST on 2010 
September 29. MBPs can be seen as intensity enhancements within the intergranular lanes. The upper- 
middle panel displays the simultaneous velocity vectors and divergence of photospheric flows. The arrows 
indicate the averaged velocity vectors determined from the local correlation tracking of G-band image 
sequences. The upper-right panel displays vorticity (s -1 , where, e.g., 0.002 s -1 represents «0.11° s -1 
in the clockwise direction) calculated from the derived horizontal photospheric flows. It is clear that 
the atmosphere in which MBPs are embedded is replete with significant horizontal flows and torsional 
motion. The lower panels, from left-to-right, display Fourier power spectra of the line-of-sight magnetic 
flux, the line-of-sight velocity and the intensity of small-scale magnetic elements (including MBPs that 
are prevalent in the upper-left panel) observed by the SOT on-board Hinode. The black circles indicate 
narrow (common) peaks that represent the observational evidence of magneto-acoustic waves in the lower 
solar atmosphere, including sausage modes, which highlight the ubiquitous nature of wave phenomena 
across all observational datasets. Images have been adapted from Morton et al. [2013] and Fujimura & 
Tsuneta [2009]. 


them to calculate the power spectrum of horizontal 
(i.e., transverse) fluctuations at the photospheric 
level. When the observational power spectrum was 
compared to a standardised Lorentzian model (see, 
e.g., the right panel of Figure 4), it was found that 
the Lorentzian power spectrum grossly underes¬ 
timated the power originating within oscillations 
>0.02 Hz (i.e., <50 s periodicity). This work sug¬ 
gests that MBP dynamics on short timescales may 
be very important in the generation of highly en¬ 
ergetic kink and/or Alfven wave phenomena. Fur¬ 
thermore, utilising both ground- and space-based 
observatories, Morton et al. [2013] presented ob¬ 
servations that revealed how kink waves can also 
be excited by the vortex motions of a strong mag¬ 
netic flux concentration in the solar photosphere 


(Figure 5). The authors detected considerable 
horizontal flows, in addition to evidence for tor¬ 
sional vorticities, in the high-resolution observa¬ 
tional ROSA data and suggested that these may 
instigate considerable wave motion in the lower at¬ 
mosphere and beyond. To test their hypothesis, 
Morton et al. [2013] computed photospheric flow 
vectors from complementary MuRAM simulations, 
which also indicated that small vortical movements 
of the photospheric plasma, with magnitudes up 
to ~0.3° s -1 , can help generate the kink waves 
observed in both high-resolution simulations and 
observations. With MBPs covering an estimated 
2.2% of quiet Sun locations [Sanchez Almeida et 
al., 2010], it does not seem inconceivable that the 
generation of kink motions in such small-scale mag- 
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netic fields may be responsible for the delivery of 
significant energy to higher atmospheric layers. 

Taking this one step further, and employing the 
long-duration observations provided by the space- 
borne Hinode/SOT instrument, Stangalini et al. 
[2013b] examined the spectral characteristics (in 
the Fourier domain) of small-scale magnetic ele¬ 
ments undergoing kink-like oscillations. Interest¬ 
ingly, the authors found that while the majority of 
the transverse oscillatory periods were in the range 
of 1—12 mHz (85—1000 s), there was no specific 
features unifying the wave phenomena originating 
within different magnetic structures. As a result, 
the authors concluded that the spectral character¬ 
istics represented a unique signature of each mag¬ 
netic element itself rather than an overarching re¬ 
lationship that defines a collective of small-scale 
magnetic structures. 

Most recently, Stangalini et al. [2014] employed 
a long time series of high-resolution photospheric 
nragnetograms to study the effects of turbulent 
convection on the excitation of kink oscillations 
in small-scale magnetic elements. Importantly, 
the authors utilised empirical mode decomposition 
techniques, which allowed them to more accurately 
analyse non-stationary time series which may be 
dominated by the horizontal displacements of mag¬ 
netic flux tubes that are continuously advected 
and dispersed by granular flows. Sub-harmonics 
of fundamental kink oscillations, with periodicities 
of 7.6 ±0.2 min, were used to verify the hypothesis 
that kink waves are induced through the buffet¬ 
ing of magnetic field lines lying at the boarders 
of photospheric convective cells. An important as¬ 
pect of this work is the potential for such kink os¬ 
cillations to be excited via non-linear interactions. 
The presence of period-doubling cascades in the 
observational results can be interpreted as a sig¬ 
nature of chaotic excitations in non-linear systems 
[e.g., Sander Sz Yorke, 2009, 2010]. This has im¬ 
portant implications for the generation of a wide 
spectrum of viable propagating periodicities in the 
lower solar atmosphere, particularly if the driving 
(i.e., buffeting) frequencies are at values that are 
not necessarily multiples or fractions of the flux 
tube’s natural frequency. Therefore, the work of 
Stangalini et al. [2014] tentatively suggests that a 
broad spectrum of wave periodicities in small-scale 
magnetic elements can be generated regardless of 
the input motions. 

3. ALFVEN WAVES 

Modern MHD simulations of the lower solar at¬ 
mosphere clearly show how torsional motions can 
easily be induced in magnetic elements in the pho¬ 
tosphere through the processes of vortical mo¬ 


tions and/or buffeting by neighbouring granules 
[Matsumoto Sz Shibata, 2010; Fedun et al., 2011b; 
Vigeesh et al., 2012]. The theoretically driven work 
of van Ballegooijen et al. [2011] and Asgari-Targhi 
Sz van Ballegooijen [2012] suggested that random 
displacements of the photospheric anchor points of 
the magnetic field lines, with velocities on the or¬ 
der of 1.5 kms -1 , would be sufficient to induce 
significant wave turbulence, thus potentially creat¬ 
ing an efficient dissipation mechanism for Alfven 
waves. Chitta et al. [2012b] provided indirect evi¬ 
dence of this effect by comparing the velocity corre¬ 
lation functions of small-scale magnetic elements in 
the photosphere. The authors detected significant 
power associated with high-frequency (>0.02 Hz; 
<50 s) horizontal motions, and suggested that 
these cases may be especially important in the cre¬ 
ation of a turbulent environment that efficiently 
promotes Alfven wave dissipation. 

Aside from the numerous theoretical studies re¬ 
lated to Alfven wave generation and dissipation, 
and as documented by Mathioudakis et al. [2013], 
observationally identifying pure torsional Alfven 
waves in the solar atmosphere has been a mon¬ 
umental struggle ever since they were postulated 
by Alfven [1942], In their most simplistic form 
(i.e., the torsional Alfven wave with azimuthal 
wave number m = 0), their incompressible na¬ 
ture provides the inability to detect them through 
typical intensity (i.e., density) measurements, and 
the azimuthal nature of the oscillation provides 
no transversal deflection of the magnetic structure 
about its central axis. Furthermore, the narrow 
nature of magnetic flux tubes in the lower solar 
atmosphere causes difficulties when attempting to 
resolve the intrinsic blue- and red-shifts associated 
with velocity measurements at opposite edges of 
the structure. However, the torsional motion of 
a magnetic element carrying an Alfven wave will 
induce a degree of non-thermal line broadening 
when observed using spectroscopic techniques [Za- 
qarashvili, 2003]. The magnitude of the broaden¬ 
ing will depend on the velocity amplitude of the 
torsional mode, which may also be compounded 
by bulk flows and/or turbulence embedded in the 
plasma. 

Employing the high-resolution spectral imag¬ 
ing capabilities of the Solar Optical Universal Po- 
larimeter [SOUP; Title et al., 1986] on the SST, 
Jess et al. [2009] were able to identify periodic non- 
thermal line broadening associated with a torsional 
Alfven wave embodied in the magnetic field lines 
anchored into a conglomeration of photospheric 
MBPs (Figure 6). Importantly, their interpreta¬ 
tion was further strengthened by the fact that the 
magnetic element did not display any periodic fluc¬ 
tuations in intensity or longitudinal/transverse ve¬ 
locity, which helped verify the absence of other 
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Figure 6. The upper-left panel displays a three-dimensional snapshot of MHD wave propagation in an 
open magnetic flux tube. The simulated image is obtained using the nonlinear magnetohydrodynamic 
SAC code, where the thin multicolor curves represent the magnetic field lines that are scaled using typical 
field strengths synonymous with MBPs (i.e., 0 — 1000 G). The upper and lower colour bars correspond 
to the magnitude of the magnetic field and the vertical velocity, 14, at the level of the photospheric 
driver, resepctively. Iso-contours of the magnetic field are displayed as solid black lines which are la¬ 
belled in the top horizontal slice taken at an atmospheric height of 1.4 Mm (or 1400 km). The bottom 
of the three-dimensional grid displays a horizontal cross-cut through the location of the torsional driver. 
The upper-right panel displays a typical expanding magnetic flux tube sandwiched between photospheric 
(broadband Ha) and chromospheric (narrowband Ha) intensity images obtained with the SST on 2007 
August 23. The observational dataset revealed that the magnetic flux tube underwent a torsional Alfvenic 
perturbation, indicated by the periodic non-thermal spectral broadening displayed in the lower panel. In 
the wavelength-versus-time plot the Ha absorption profile line width, calculated using spectral imaging 
techniques, is observed to oscillate with a periodicity ~420 s, with consecutive peaks indicated by the 
white arrows. As can be seen in both upper panels, the Alfvenic displacements are torsional motions 
that remain perpendicular to both the direction of propagation and the magnetic fields outlining constant 
magnetic surfaces. Images adapted from Jess et al. [2009] and Fedun et al. [2011b]. 


magneto-acoustic modes. In follow-up work, Math- 
ioudakis et al. [2013] revealed the independent and 
opposite Doppler shifts associated with the oppos¬ 
ing edges of the magnetic element, thus reinforc¬ 
ing the interpretation that the observational sig¬ 
natures represented a torsional Alfven wave. In 
this work, it was found that the amplitude of the 
non-thermal broadening was ~0.05 A, equating to 
a velocity amplitude of ~2.5 kms” 1 . When com¬ 
bined with a local Alfven speed on the order of 


22 kms” 1 and a plasma density of approximately 
1 x 10” 9 gem” 3 , the resulting wave energy was 
estimated as 150 000 Wm 2 . While this is a 
vast quantity of available energy, far in excess of 
the threshold required to sustain localised chromo¬ 
spheric and coronal heating, the true importance 
can only be ascertained once the filling factor of 
such waves is accurately known (see, e.g., Chap¬ 
ter 24). However, as documented by Goossens et 
al. [2011] and Mathioudakis et al. [2013], Alfven 
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waves are naturally difficult to dissipate unless 
they are able to find an alternative mechanism to 
promote efficient energy dissipation. Such mecha¬ 
nisms include phase mixing and resonant absorp¬ 
tion [e.g., Goossens & De Groof, 2001], whereby 
non-uniformities in the magnetic field configura¬ 
tions results in the coupling between neighbour¬ 
ing magnetic iso-surfaces, thus promoting a signif¬ 
icantly steep gradient to allow efficient energy dissi¬ 
pation. Other possibilities are the mode conversion 
of Alfven waves into magneto-acoustic modes that 
can propagate obliquely to the magnetic field lines, 
thus promoting efficient energy loss [e.g., Parker, 
1991; Nakariakov et ah, 1997], or turbulent mix¬ 
ing as a result of high frequency fluctuations [e.g., 
van Ballegooijen et al., 2011; Asgari-Targhi & van 
Ballegooijen, 2012], 

Utilising the same dataset as presented by Jess et 
al. [2009], Fedun et al. [2011b] examined the spatial 
structuring of the observed torsional Alfven fre¬ 
quencies and related these to the outputs of nonlin¬ 
ear three-dimensional magnetohydrodynamic nu¬ 
merical simulations from the Sheffield Advanced 
Code [SAC; Shelyag et al., 2008]. The authors im¬ 
plemented a vortex driver at the base of the simu¬ 
lated domain (see the left-hand image of Figure 6) 
and revealed how magnetic flux tubes can act as a 
spatial frequency filter for torsional Alfven waves. 
Importantly, the authors found that his form of 
frequency filtering is strongly dependent on the 
structure and geometry of the magnetic field it¬ 
self. This implies that the observed spatial wave 
power and oscillatory frequencies can be a function 
of the underlying MBP, possibly allowing magnetic 
fields to be mapped as a function of atmospheric 
height solely employing such seismology techniques 
[Fedun et al., 2011b], 

Even with higher sensitivity equipment now be¬ 
coming commonplace on a majority of ground- 
based telescopes [e.g., the IBIS and CRISP imag¬ 
ing spectropolarimeters; Cavallini, 2006; Scharmer 
et al., 2008], there is still a significant lack of sub¬ 
sequent Alfven wave detections in the lower so¬ 
lar atmosphere. This is in stark contrast to mod¬ 
ern numerical simulations [e.g., the MuRAM code; 
Vogler et al., 2005] which indicate widespread tor¬ 
sional motions in synthesised photospheric filter- 
grams [Shelyag et al., 2013]. Instead, many have 
turned their attention to the chromosphere where 
the interaction, mode conversion and dissipation 
of Alfven waves may have more identifiable sig¬ 
natures, particularly in regions where the plasma 
/3 changes abruptly, or self-induced turbulence re¬ 
sults in rapid localised dissipation [van Ballegooi¬ 
jen et al., 2011; Asgari-Targhi & van Ballegooijen, 
2012 ], 


4. CONCLUSIONS & FUTURE DIRECTIONS 

In the near future it is expected that high 
sensitivity 2D spectropolarimeters (e.g., CRISP 
and IBIS) will be employed simultaneously along¬ 
side high-cadence imagers (e.g., Hinode/SOT and 
ROSA) to obtain multiwavelength time series at 
the highest spatial, temporal and spectral reso¬ 
lutions currently achievable (see, e.g., Figure 7). 
High precision measurements will allow the char¬ 
acterisation of MHD waves themselves manifesting 
in small-scale magnetic elements that are at the 
limits of current telescope resolving power. Im¬ 
portant oscillatory parameters, such as the prop¬ 
agation speeds, amplitudes and phase relation¬ 
ships will allow MHD wave phenomena to be doc¬ 
umented with unprecedented accuracy, including 
the establishment of evidence to verify the presence 
of standing/propagating, fast/slow, trapped/leaky 
and surface/body oscillatory modes. Furthermore, 
the multiwavelength nature of the data will also 
enable the detected MHD waves to be tracked 
through the solar atmosphere as they journey from 
the photosphere, through the tenuous chromo¬ 
sphere, and into the super-heated corona. Impor¬ 
tantly, coverage of the waves as they propagate 
through the different atmospheric regions where 
the plasma f3 changes from magnetically domi¬ 
nated to plasma pressurised regimes will provide 
valuable insight on aspects of mode coupling and 
wave dissipation. Many previous examples have in¬ 
dicated that regions where /3 = 1, often in locations 
sandwiched between the photosphere and chro¬ 
mosphere, provides opportune atmospheric condi¬ 
tions to promote efficient oscillatory mode conver¬ 
sion, thus allowing naturally difficult-to-dissipate 
waves (e.g., Alfven waves) to convert into more 
readily dissipated compressible modes [e.g., Ulm- 
schneider et al., 1991; Kalkofen, 1997; Hasan et al., 
2003; McAteer et al., 2003; Bloomfield et al., 2006; 
Goossens et al., 2006; Jess et al., 2012b, to name 
but a few]. 

In addition to maximising the scientific return of 
the current fleet of ground- and space-based solar 
telescopes, the solar community eagerly awaits the 
arrival of the first next-generation high-resolution 
facilities, including the 2 m National Large Solar 
Telescope [NLST; Hasan et al., 2010] in Ladakh, 
India, and the 4 m Daniel K. Inouye Solar Tele¬ 
scope [DKIST, formerly the Advanced Technology 
Solar Telescope, ATST; Keil et al., 2003; Rim- 
mele et al., 2010] atop the Haleakala volcano on 
the Pacific island of Maui. These pioneering facil¬ 
ities are due to receive first light towards the lat¬ 
ter stages of this decade, and by utilising dramat¬ 
ically increased aperture sizes (4 m in the case of 
DKIST), photospheric structures down to ~20 km 
in size will be able to be detected, tracked and stud- 
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Figure 7. Simultaneous images of a solar active region, captured using both existing and next-generation 
imaging detectors on the DST during 2014 August 24. The left and middle images are G-band (photo¬ 
sphere) and Ca II K (upper-photosphere/lower-chromosphere) snapshots, respectively, acquired using the 
electron-multiplying CCD cameras of the ROSA multi-camera system. A true-size representation of the 
Earth is depicted in the lower-right section of each image to provide a sense of scale. The image on the 
right reveals the solar chromosphere, captured through a 0.25 A Ha-core filter, and employs a 4.2 MP 
Zyla CMOS sensor from Andor Technology, which allows image sequences to be obtained with 15 ms 
exposure times and frame rates exceeding 60 s -1 . Furthermore, the large pixel formats of modern CCD 
and CMOS detectors allow fields-of-view in excess of 200” x 200” to be sampled at the diffraction limit, 
and when combined with numerous detectors each sampling discreet wavelengths, provides a seamless 
view through the tenuous atmospheric layers. Combining such high-sensitivity detectors with current 
and future telescope facilities will open the door for greater scientific understanding through drastically 
improved number statistics and larger fields-of-view. Images courtesy of D.B. Jess. 


ied in unprecedented detail. In the longer term 
is the launch of the Solar-C payload (expected 
launch date around 2019) with a UV/visible/IR 
telescope approaching 1.4 m in diameter [Suematsu 
et ah, 2014], The seeing free data from a high- 
resolution space-based facility will enable near- 
continuous observations of magnetic elements to 
be obtained without the need for extensive post¬ 
processing algorithms. Even more long-term is the 
4 m European Solar Telescope [EST; Collados et 
al., 2010] and the 5-8 m Chinese Giant Solar Tele¬ 
scope [CGST; Liu et ah, 2014], which are expected 
to begin construction following 2020, and allow 
structures as small as ~10 km in size to be ob¬ 
served for the first time. 

The desire to improve not only the spatial res¬ 
olution of solar telescopes, but also the tempo¬ 
ral and spectral resolutions through the develop¬ 
ment of new imaging and spectropolarimetric in¬ 
struments, will enable key outstanding questions 
related to lower atmospheric understanding to be 
firmly addressed. Such questions include: 

• Do all small-scale magnetic elements carry 
energy in the form of MHD waves? There is 
an abundance of evidence that demonstrates the 
suitability of narrow magnetic strands as effi¬ 
cient wave/energy conduits. For example, Jess et 
ah [2012c] found that ~73% of MBP structures 
demonstrated magneto-acoustic wave signatures. 
However, does this mean that ~27% of these fea¬ 


tures contain no oscillatory phenomena? Or do 
waves exist, but under different guises (e.g., the 
more difficult to detect Alfven modes)? 

• What is the dominant mode of oscillation 
present in small-scale photospheric magnetic fea¬ 
tures? The majority of research to date has 
focussed on the magneto-acoustic signatures de¬ 
tectable through intensity (i.e., density) fluctu¬ 
ations. However, with improved observations 
and detection algorithms [e.g., Stangalini et ah, 
2013a, b], it has become clear that other wave 
modes may be present in small-scale magnetic ele¬ 
ments alone or alongside their compressible coun¬ 
terparts. Are these additional modes (e.g., kink 
and Alfven waves) superimposed on top of the 
seemingly ubiquitous compressive oscillations? Or 
are they the dominant mode of oscillation in cer¬ 
tain magnetic features? And if so, which features? 

• How do the (5 = 1 layers contribute to both 
the waveforms visible in the outer solar atmosphere 
and the rate of localised heating through wave dis¬ 
sipation? It is well accepted that regions of the 
Sun’s atmosphere where (5 = 1 (i.e., the mag¬ 
netic pressure is equal to the gas pressure) are op¬ 
timal for oscillatory mode conversion [e.g., Ulm- 
schneider et ah, 1991; Kalkofen, 1997; Hasan et 
ah, 2003]. McAteer et ah [2003] and Bloomfield et 
ah [2006] have demonstrated the coupling between 
kink waves and longitudinal oscillations, while Jess 
et ah [2012b] have documented the conversion be¬ 
tween longitudinal waves and kink/sausage modes. 
Regions where (5 = 1 are also predicted to have 
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important consequences for Alfven wave conver¬ 
sion through processes such as resonant absorption 
[e.g., Goossens et ah, 2006]. Can such atmospheric 
layers convert observed (and possibly unobserved!) 
waves into other modes that can readily dissipate 
their energy into the surrounding plasma? Do the 
preferred conversion mechanisms substantiate the 
abundance of oscillatory motion observed in the 
outer regions of the solar atmosphere? 

• Can we find, evidence for the existance of 
magneto-acoustic wave modes with azimuthal wave 
numbers exceeding 1? The observational work to 
date suggests an abundance of sausage (m = 0) 
and kink (rri =1) mode oscillations in the lower 
solar atmosphere. However, theoretical predictions 
also suggest that higher order wave numbers (i.e., 
fluting modes with m > 2) should also be preva¬ 
lent in magnetic flux tubes [for a recent review see 
Goossens et ah, 2011, and the references therein]. 
Why do we currently have no evidence for such os¬ 
cillations? Is it a result of relatively poor telescope 
resolution and/or instrument sensitivity? Or are 
the amplitudes of these oscillatory modes so small 
that they become impossible to disentangle from 
other superimposed waves? 

• What role do downwardly propagating MHD 
waves play in the structuring, evolution, dynamics 
and energy balance of magnetic features in the solar 
photosphere? Attempting to find the elusive solu¬ 
tion to the coronal heating problem has resulted in 
an overwhelming interest in upwardly propagating 
MHD waves. However, observations continue to 
identify significant oscillations propagating down¬ 
wards through the solar atmosphere [e.g., Gupta et 
al., 2013]. Jess et al. [2012c] estimated that approx¬ 
imately a quarter of all detected magneto-acoustic 
waves in MBPs were downwardly propagating. Are 
these propagating oscillations generated at higher 
atmospheric heights? Are they the counterpart 
of upwardly propagating oscillations generated at 
the opposite foot point of a closed loop? Or are 
they the result of reflection as upwardly propagat¬ 
ing waves encounter the severe density discontinu¬ 
ities intrinsic to chromospheric and transition re¬ 
gion layers? Furthermore, how do these downward 
motions contribute to the mass/energy flow asso¬ 
ciated with magnetic elements in general? 

The above questions are not intended to high¬ 
light all of the outstanding problems related to 
MHD wave phenomena in the lower solar atmo¬ 
sphere, nor are they listed in order of importance. 
Instead, they are simply listed as those which we 
feel have an overarching central importance when 
attempting to address the long-standing issues of 
wave generation, propagation, energy transfer and 
dissipation throughout the photosphere and be¬ 
yond. Furthermore, it seems unlikely that the 
above questions can be unequivocally answered re¬ 


lying solely on observational approaches. A com¬ 
bination of theoretical, analytical and numerical 
modelling techniques will be required to help ex¬ 
tract and interpret the wealth of MHD wave modes 
that exist in the lower solar atmosphere. Thus, 
the answers to the key science questions outlined 
above will only arise through the novel use of high- 
resolution (spatial, temporal and spectral) photo- 
spheric datasets alongside the rapid development 
of our theoretical understanding. 
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